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1.  INTRODUCTION 


The  laser  ignition  concept  for  initiating  large-calibCT  ammunition  was  originally 
conceived  and  tested  at  the  U,S.  Ballistic  Research  Laboratory*  at  Aberdeen  Proving 
Ground,  MD  21005-5066  (Robitalle  1964;  Barrows  et  al.  1993).  The  ability  to  ignite 
propellant  beds  by  using  laser  radiation  as  the  ignition  source  should  eliminate  the  need  for 
primers  and  simplify  the  ignition  train  of  ammunition  as  well  as  improve  the  safety  of  the 
firing  procedure. 

However,  if  the  laser  radiation  cannot  reach  the  igniter  material  that  subsequently 
ignites  the  propellant,  the  propellant  cannot  ignite.  Important  to  all  schemes  that  transmit 
laser  radiation  into  the  breech  of  a  gun  is  an  optical  window.  This  window  must  be  maHp. 
of  a  material  that  transmits  the  laser  radiation  without  deletaious  effects  to  the  material  and 
must  withstand  the  high-tercperature,  high-pressxire,  and  chemically  hostile  environment  of 
the  breech  during  the  ballistic  cycle.  In  addition,  the  window  material  and/or  design  must 
eschew  particles  created  during  the  ballistic  cycle  that  would  otherwise  be  deposited  on  the 
window.  If  particles  remain  on  the  window  surface,  the  laser  radiation  transmission  will 
be  reduced  for  the  subsequent  firing  (see  Figure  1).  At  some  point,  if  the  window  were  not 
sufficiently  cleaned,  it  would  no  longer  transmit,  and  further  firing  would  be  impossible. 
Also,  if  thermal  shock  fix)m  contact  with  the  hot  gases  were  to  occur,  the  window  would 
crack  and  the  pressure  seal  would  be  lost 


On  30  September  1992,  the  U.S.  Army  Ballistic  Research  Laboratory  (BRL)  was 
deactivated  and  subsequently  became  part  of  tiie  U.S.  Army  Research  Laboratory  (ARL) 
on  1  October  1992. 
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Igniter  Pad 


Figure  1.  Effect  of  dirtv  window  on  laser  ignition. 


A  cleaning  procedure  that  required  swabbing  the  window  after  each  firing  would 
not  be  acceptable  for  many  reasons,  the  least  of  which  is  the  requirement  for  rapid  firing, 
aeaning  with  a  bmsh  mounted  on  the  breech  should  meet  with  only  partial  success.  The 
brush  would  rapidly  become  contaminated  and  in  a  steady-state  condition  would  leave 
particles  on  the  window  as  it  brushes  others  away.  Therefore,  a  breech  brash  would 
require  Sequent  replacement  in  order  to  avoid  the  steady-state  condition. 


A  cleaning  procedure  that  would  work  with  only  a  subset  of  the  ammunition  (those 
that  require  a  stub  base)  would  be  a  double  window.  If  a  window  were  mounted  in  the 
stub  base  as  well  as  in  the  breech,  the  stub  base  window  should  provide  an  optically  clean 
screen  for  the  breech  window  that  would  be  thrown  away  with  every  round.  Such  a 
procedure  would  dramatically  increase  the  cost  of  ammunition  since  acceptable  windows 
are  not  inexpensive.  This  paper  presents  a  change  in  paradigm  fiom  the  aforementioned 
methods.  A  technique  that  uses  the  gases  fi'om  propellant  combustion  to  keep  the  window 
clean  is  discussed.  Hydrodynamic  forces  will  always  be  present  in  a  gun  as  long  as 
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propellant  and  hot,  high-pressure  gases  are  used  to  propel  the  projectile.  If  a  cavity  is 
placed  in  the  spindle,  these  forces  can  be  utilized  to  keep  the  window  clean. 


2.  CALCULATION  OF  PRESSURE  IN  THE  CAVITY 


Figure  2  schematically  represents  the  flow  system  composed  of  the  cavity  in  the 
spindle  and  the  gun  chamber.  The  two  flow  volumes  are  connected  by  two  orifices  of 
equal  diameter.  During  the  early  phase  of  the  intmor  ballistic  cycle,  the  pressure  in  the  gun 
chamber  rises  following  the  combustion  of  the  propelling  charge.  A  flow  of  propellant 
gases  to  the  cavity  through  the  orifice  then  occurs.  At  some  point  in  the  ballistic  cycle  after 
the  gun  chamber  pressure  has  reached  its  peak  value,  the  pressure  in  the  cavity  becomes 
higher  than  that  in  the  gun  chamber.  The  flow  reverses  its  direction  until  the  end  of  the 
interior  ballistic  cycle. 


Figure  2.  Proposed  flow  system  for  hvdrodvnamic  cleaning  of  window. 


Since  the  gas  flow  through  the  orifice  is  a  flow  with  high  Reynolds  numbers,  it  is 
appropriate  to  consider  the  flow  as  a  one-dimensional,  inviscid  flow  system  for  the  primary 
interest  of  predicting  the  pressure  rise  in  the  cavity.  The  following  sections  contain  the 
governing  equations  for  the  flow  system. 
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(1) 


2.1  Equation  of  State.  The  Nobel- Abel  equation  of  state  is  adopted  and  written  as: 

p  Ru  Tc  W 

Vc  m  -  b  W  ’ 


where  Pc,  Tc,  W,  and  Vc  are  the  pressure,  tenapwature,  total  mass,  and  volume  of  the 
cavity,  respectively.  Furthermore,  m ,  and  Ru  are  the  molecular  weight,  covolume,  and 
universal  gas  constant  of  tiie  propellant  gases,  respectively. 


2.2  Continuity  Equations.  It  is  noted  that  there  are  two  modes,  choked  and 
nonchoked,  of  flow  through  the  orifice  during  the  interior  ballistic  cycle.  A  choked  flow 
occurs  when 


P*  r  2  ,m-l) 
Fo  “  V 


where  P*  =  pressure  of  expanding  fluid  at  which  a  choked  flow  occurs,  Po ,  To  equal  the 
pressure  and  temperature  in  gun  chamber,  and  k  =  ratio  of  specific  heats.  For  propellant 
gases,  k  =  1.243,  and  the  pressure  ratio  in  equation  (2)  is  calculated  to  be  0.556.  In  the 
equation,  Po  can  be  obtained  fiom  calculations  using  XKTC  code  or  others  (Gough  1986). 


2.2. 1  Flow  From  the  Gun  Chamber  to  the  Cavity.  When  ^  <  0.556  (choked 
flow):  The  flow  rate  through  the  orifice  is  (Shapiro  1953) 

where  Cj  =  discharge  coefficient  of  orifice,  A  =  cross-sectional  area  of  orifice,  and  g  = 
gravity. 


When 


To 


>  0.556  (nonchoked  flow): 


The  flow  rate  is 


(4) 


where  M  =  Mach  number  given  as 


p 

2.2.2  Flow  from  the  Cavity  to  the  Gun  Chamber.  When  <  0.556  (choked 
flow):  The  flow  rate  is 

rA.#+l)/»-l) 

p 

"When  ^  >  0.556  (nonchoked  flow);  The  flow  rate  is 


with 


(7) 


(8) 


2.3  Energy  Equations. 

2.3.1  Flow  From  the  Gun  Chamber  to  the  Cavity.  For  either  a  choked  or  nonchoked 
flow,  the  rate  of  temperature  change  in  the  cavity  is 

(O 

[(1”  To  ■  T’c]  ,  (9) 

where  ^  =  heat  loss  to  the  walls  of  the  cavity,  including  the  orifice  area;  Cp  =  constant 
pressure  specific  heat;  and  Cy  =  constant  volume  specific  heat. 

2.3.2  Flow  From  the  Cavity  to  the  Gun  Chamber. 

(oT  r 

=  wCv  ^  "^v)  •  (10) 
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3. 


EXPERIMENTAL 


3.1  Water  Jet  Simulator.  Initial  visualizing  of  the  diverted  jet  concept  was 
performed  using  water.  Figure  3  shows  a  water  tank  with  provisions  for  pressurized  water 
addition  as  well  as  a  removable  bottom  plate  that  contained  the  jet  orifices.  Plates  were 
tried  with  angles  of  30°,  45°,  and  60°  (see  Figine  2)  between  the  jets.  The  larger  angles 
required  a  two-stage  structure  for  the  diverting  jet  orifice.  Otherwise,  at  large  angles  the 
entrance  to  the  diverting  jet  would  be  too  far  from  the  direct  jet  Therefore,  a  large  entrance 
was  drilled  part  way  through  the  plate  and  was  parallel  to  the  direct  jet  The  diverting  jet 
was  then  drilled  so  as  to  intersect  the  entrance  hole  and  the  exit  of  the  direct  jet 


Figure  3.  Schematic  of  water  iet  visualization  system. 


(see  Figure  2).The  tank  system  was  supported  by  four  acrylic  tubes  over  a  drainage 
system.  Operation  of  the  system  was  then  videotaped. 

3.2  Flamespread  Chamber.  For  subscale  tests,  the  flamespread  chamber 
(Figure  4)  was  used  (Kooker,  Chang,  and  Howard  1992;  Kooker,  Howard,  and  Chang 
1994).  The  flamespread  chamber  consisted  of  a  transparent  acrylic  tube  (interior  diameter 
of  76  mm  with  an  axial  dimension  of  305  mm)  contained  in  a  steel  confinement  casing. 
Ports  were  machined  in  the  steel  casing  and  the  acrylic  tube  for  pressure  transducers. 
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Figure  4.  Cross-sectional  view  of  flamesuread  simulator  with  cavity  and  igniter  chambers. 

For  these  experiments,  the  acrylic  mbe  was  split  into  two  regions  (igniter  and 
breech  window  regions).  The  igniter  region  of  the  tube  was  closest  to  the  right  side  of  the 
simulator  and  contained  a  cloth  basepad  (containing  black  powder,  ball  powder,  or  clean 
burning  [CBI]  material).  The  basepad  was  ignited  by  a  low-shock  igniter  located  on  the 
outside  of  the  basepad  facing  the  breech  region  so  that  laser  ignition  could  be  modeled  (see 
Figure  4).  A  pressure  transducer  was  located  in  each  chamber.  Other  pressure  transducers 
were  located  in  various  parts  of  the  chambCT,  simulating  the  laser  breech  window  area.  The 
low-shock  igniter  was  electrically  connected  to  the  firing  line  via  a  high-pressure  electrical 
feedthrough  in  the  top  of  the  simulator.  The  chamber  was  also  placed  on  the  hori2X)ntal  to 
approximate  typical  tank  or  artillery  use. 


The  window  region  was  experimentally  modelled  by  a  4.25-mm  O.D.  aluminum 
plug  (a  threaded  right-circular  cylinder)  that  was  highly  polished  on  one  end  to  form  a 
wimess  plate  and  slotted  on  the  other  end  so  that  it  could  be  screwed  out  of  the  hollow  tube 
for  viewing  (the  hollow  tube  was  made  of  aluminum  and  was  internally  threaded  at  the 
window  region  so  that  it  could  tightly  hold  the  witness  plate  in  place).  The  plug  was 
chosen  so  that  a  leak-proof  seal  could  be  maintained  in  the  window  region  and  yet  allow 
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easy  access  to  the  witness  plate  for  evidence  of  residues  from  the  igniter  material 
combustion. 

3.3  Full-Scale  Simulator.  The  full-scale  simulator  was  a  165-mm  interior  diameter 
by  1,000-mm-long  acrylic  tube  with  aluminum  end  plates  held  in  place  by  four  threaded 
rods  (see  Figure  5).  The  igniter  chamber  was  scaled  up  from  the  subscale  simulator 
chamber.  The  conical  cavity  chamber  was  cast  in  the  bottom  of  the  acrylic  tube  with  an 
aluminum  cover  plate  that  contained  the  jet  orifices  and  isolated  the  cavity  from  the  main 
chamber  that  contained  the  basepad.  The  casting  material  was  a  pourable  polystyrene- 
polyester  mix  that  was  transparent  upon  sohdification.  Operation  of  the  full-scale  simulator 
was  identical  to  that  of  the  subscale  simulator. 


Figures.  Schematic  of  full-scale  simulator. 


4 .  RESULTS  AND  DISCUSSION 

4.1  Calculations.  The  governing  equations  described  previously  can  be  integrated 
numerically  using  a  conventional  integration  method  (see  Appendix).  The  input  data  for  the 
computer  program  are: 
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Discharge  coefficient  of  orifice  Cd  =  0.62 
Cross-sectional  area  of  the  orifice  A  =  0.0000565 
Ratio  of  specific  heats,  k  =  1.243 
Volume  of  cavity,  Vq  =  0.0000872  and  0.000127  m3 
Specific  covolume,  b  =  0.00103  m3/kg 
Molecular  weight,  m  =  23.4  kg/kg-mole 
Constant  volume  specific  heat,  Cy  =  456  J/K-kg 
Constant  pressure  specific  heat,  Cp  =  566  J/K-kg 
Gravitational  acceleration,  g  =  9.81  m/s^ 

Conversion  factor,  J  =  778 
Universal  gas  constant,  Ru  =  83 14  J/K  kg-mole 
Initial  pressure  in  the  cavity,  Pco  =  0.101  MPa 
Initial  temperature  in  the  cavity,  Tco  =  295  K 
Heat  loss  to  the  cavity  walls,  (p  =  0A  and  0.2  . 

Note  that  the  cross-sectional  area.  A,  given  here  is  the  combined  cross-sectional  area  of  the 
two  orifices  in  the  window  design,  each  with  a  diameter  of  4.25  mm  (0.165  inch).  This 
diameter  is  sufficient  for  a  laser  beam  to  pass  through.  It  is  possible  to  open  a  cavity  inside 
the  spindle  with  a  volume,  Vc,  ranging  from  0.0000872  to  0.000127  m3  without  causing  a 
stmctural  failure  of  the  spindle  upon  high-pressure  loadings  during  the  ballistic  cycle.  As 
for  heat  loss  to  the  cavity  walls,  (p  is  difficult  to  accurately  determine  without  conducting  a 
complex  three-dimensional  flow  analysis.  However,  from  experience  gained  from  closed 
bomb  testing  of  propellants,  it  is  reasonable  to  assume  the  heat  loss  to  be  in  the  range  of  10 
to  20%  of  the  total  energy  entering  the  two  orifices.  In  a  continuous  firing  series,  the  heat 
loss  will  decrease  accordingly. 


A  smaller  Vc  will  result  in  a  higher  pressure  rise  in  the  cavity.  Figure  6  presents 
calculated  pressure  rises  in  the  cavity  with  Vc  =  0.000127  m3  corresponding  to  two 
assumptions  of  heat  loss.  During  the  interior  ballistic  cycle,  the  maximum  pressure 
differentials  between  the  gun  chamber  and  the  cavity,  AP  =  Po-  Pc,  are  180  MPa  and 
2(X)  MPa  for  (p  =0.10  and  <p  =  0.20,  respectively.  Based  on  the  larger  value  of  these  two 
pressure  differentials,  the  shear  stress,  S,  across  the  solid  between  the  gun  chamber  and  the 
cavity  is  calculated  by  Equation  11. 


AP  Ai 
~  A2 


(11) 


where  A  i  is  the  surface  area  on  which  the  pressure,  AP,  applies  and  A2  is  the  cross- 
sectional  area  across  the  solid.  The  values  of  Aj  and  A2  for  the  proposed  cavity  are 
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0.00310  and  0.00501  respectively.  The  resulting  maximum  shear  stress  on  the 
interface  between  the  gun  chamber  and  the  cavity  chamber  is  124  MPa. 


■ -  Gun  Chamber 

inmamiii  10%  Heat  Loss 
20%  Heat  Loss 

;xxxxxx;>  Ap  10% 

-  Ap20% 


Figure  6.  Pressure  and  Ap  curves  of  gun  chamber  and  cavity  chamber  with  heat  losses  of 
10%  and  20%. 

4.2  Water  Jets.  The  water  tank  was  connected  to  the  domestic  water  supply  at  a 
pressure  of  approximately  0.3  MPa  (40  psig).  The  tank  was  vented  and  the  water  flow 
from  the  exit  orifices  observed.  The  flow  was  diverted  at  an  angle  from  the  direct  jet  orifice 
(see  also  Figures  2  and  3).  The  diverting  angle  was  slightly  less  than  half  the  angle  of  the 
diverting  jet.  This  result  was  expected.  At  a  diverting  jet  angle  of  60°,  the  diverted  distance 
was  greatest  for  the  angles  tested.  The  water  jet  appeared  to  completely  miss  the  region 
where  the  window  would  be  placed  (see  the  object  under  the  jets  in  Figiure  3).  The  angle  of 
60°  was  the  largest  angle  possible  for  the  machining  techniques  utilized.  Therefore,  this 
angle  was  chosen  for  further  study  in  the  simulators. 

4.3  Flamespread  Chamber  (Subscale  Simulator).  The  first  tests  with  the  subscale 
simulator  used  a  simple  chamber  (a  right-circular  cylinder)  for  the  cavity  (see 

Figure  4).  A  baseline  test  with  Class  3  black  powder  was  run  with  only  a  direct  jet  (e.  g., 
the  diverting  jet  was  not  completely  drilled)  and  without  the  cavity  tube  shown  in  Figure  4. 


10 


The  black  powder  residue  particles  were  captured  on  a  mylar  sheet  placed  over  the  window 
orifice.  Comparison  of  the  before  and  after  photographs  in  Figure  7  of  the  orifice  entrances 
showed  that  black  powder  residue  in  the  gun  chamber  portion  of  the  simulator  is  indeed 
very  dirty.  As  shown  in  Figure  8,  the  number  of  particles  that  passed  through  the  direct  jet 
was  quite  large.  It  is  doubtful  that  a  laser  beam  penetrating  this  residue  layer  (without  first 
a  cleaning  procedure)  could  successfully  ignite  an  igniter  pad. 


In  addition  to  not  directly  hitting  the  window  region  by  the  opening  of  the  diverting 
jet  for  the  remainder  of  the  tests,  the  ignite’  gases  were  prevented  from  directly  contacting 
the  window  witness  plate  by  placing  a  hollow  tube  of  4.25-mm  LD.  tube  that  extended 
from  the  window  witness  plate  to  nearly  the  exit  of  the  orifices  that  form  the  jets  (see 
Figures  2  and  4).  It  was  diought  that  the  initial  air  in  the  tube  would  act  as  a  "buffer"  that 
would  compress  as  ignite  gases  entered  the  cavity  chamber,  thereby  preventing  a  direct 
path  to  the  window  for  both  particles  and  hot  gases.  These  tests  also  used  Class  3  black 
powder  in  the  basepad. 

This  pardcular  technique  drastically  reduced  the  particles  that  attached  to  the 
window  region.  Howevar,  protection  of  the  window  was  not  particularly  successful.  A 
slight  film  and  a  large  number  of  particles  still  covered  the  window  wimess  plate  (see 
Figure  9b).  Fortunately,  it  was  noted  that  the  majority  of  particles  in  the  diverted  gas  flow 
hit  the  chamber  wall  and  then  appeared  to  flow  back  toward  the  front  of  the  cavity.  It  was 
posited  that  this  recirculating  flow  carried  die  particles  that  finally  attached  to  the  window. 
Borrowing  technology  commercially  available  for  separating  particles  from  gas  flows,  the 
cavity  chamber  was  changed  from  a  right-circular  cylinder  to  a  tapered  cone  with  the 
minimum  diameter  in  the  window  region.  This  shape  formed  a  hydrocyclone  (see 
Figure  10a)  that  would  separate  particles  from  a  tangential  gas  flow.  While  the  diverted 
flow  into  the  cavity  chamber  is  not  purely  tangential,  the  flow  becomes  largely  tangential 
tq)on  impacting  the  curved  wall  of  the  cone.  The  tangential  velocity  component  increases  if 
the  diverting  jets  are  offset  from  the  center  of  the  cone.  Fortunately,  such  is  the  case  for  the 
breech  spindle  currently  in  use  for  155-mm  caimon-the  primer  hole  is  offset  by  14  mm  to 
accommodate  such  features  as  die  Swedish  notch. 
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(b) 

Figure  10.  Hvdrocvclone  chamber  (a')  and  window  witness  plate  fbV 

This  technique  reduced  to  single  digits  the  number  of  particles  that  finally  found  the 
window  (see  Figure  10b).  When  clean-burning  igniter  (CBI)  material  instead  of  black 
powder  was  used,  particles  were  not  noticed  on  the  window  witness  plate.  It  was  decided 
that  black  powder  represented  the  worst  case.  An  experiment  was  then  tried  with  ten 
consecutive  shots  of  black  powder  before  examining  the  witness  plate.  The  wimess  plate 
had  few  particles  and  was  only  slightly  filmed.  A  larger  number  of  shots  would  have  been 
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needed  before  laser  light  would  begin  to  be  seriously  attenuated.  Fortunately,  at  low  rates 
of  loading  of  the  window  the  laser  can  bum  off  some  of  the  contaminates  as  it  ignites 
successive  shots,  thereby  increasing  the  total  number  of  shots  before  cleaning  of  the 
window. 


4.4  Full-Scale  Simulator.  The  full-scale  simulator  was  made  to  resemble  the  end 
of  the  breech  chamber  and  the  spindle  face  in  a  155-mm  cannon.  The  simulator  utilizing 
the  conical  cavity  chamber  showed  no  deviation  from  the  results  of  the  subscale  simulator. 
As  evidenced  by  residue  particles  that  attached  themselves  to  the  wall  of  the  conical 
chamber,  the  gas  flow  from  the  diverting  jet  avoided  the  window  region  and  appeared  to 
move  in  a  circular  path  in  the  cavity  chamber.  The  window  region  remained  essentially 
clean  (similar  to  the  witness  plate  in  Figure  10b)  while  using  black  powder  as  the  igniter 
material. 


5.  SUMMARY 


Both  liquid-  and  gas-phase  jets  exiting  fiom  a  pressurized  chamber  into  a  second 
gas-filled  chamber  of  lower  pressure  can  be  diverted  from  the  line-of-sight  through  the  jet 
orifices.  This  concq)t  was  used  in  a  155-mm  cannon  simulator  to  reduce  the  dirt  and 
thermal  shock  that  can  affect  the  window  used  for  laser  ignition.  Gases  emanating  fi-om  an 
igniting  basepad  that  would  otherwise  impinge  upon  and  obscure  the  window  were 
diverted.  The  gases  were  further  processed  by  centrifugal  force  within  a  conical  cavity  next 
to  the  window  to  remove  solid  particles  from  the  gas  flow  that  could  otherwise  attach  to  the 
window.  An  air  buffer  zone  just  above  the  window  was  also  provided.  These  two 
techniques  kept  hot,  dirty  gases  from  the  window.  They  should  be  effective  in  protecting 
the  window  fi-om  thermal  shock  and  fi-om  becoming  dirty,  which  would  block  the  laser 
beam  required  to  ignite  the  basepad. 
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APPENDIX: 

FORTRAN  PROGRAM  CODE 


17 


There  are  three  files: 

sourceT  =  source  code  listing 
data.inc  =  include  file 
constants.inc  =  include  file 

"Include"  files  are  pieces  of  code  which  get  inserted  into  the  program  at  compilation, 
Constants  and  common  block  information  were  placed  into  these  files. 

To  compile  the  code,  SLATEC  libraries  must  be  installed  on  the  system. 

Til  -o  chamber  source.f  /usr/lib/slatec/lib/slatec.a 

The  SLATEC  (acronym  for  Sandia,  Los  Alamos,  Air  Force  Weapons  Laboratory  Technical 
Exchange  Committee)  mathematic  libraries  perform  a  wide  variety  of  functions  such  as 
solving  linear  systems  of  equations,  ordinary  and  partial  differential  equations,  spline 
interpolation,  etc.  This  FORTRAN  library  is  in  the  public  domain  and  can  be  obtained 
from  the  following  organization: 

Energy  Science  and  Technology  Software  Center 

P.O.  Box  1020 

Oak  Ridge,  TN  37831 

Telephone  615-576-2606 

E-mail:  estsc%a1.adoitis.mrouter@zeus.osti.gov 

The  following  codes  were  compiled  on  an  Silicon  Graphics  Indigo.  Silicon  Graphics  Inc. 
freely  distributes  the  SLATEC  libraries  for  their  workstations  and  file  servers. 


Filename:  source.f 


PROGRAM  CHAMBER 


C - 

C  State  Variables: 

C 

C  PC  =  chamber  pressure  (psi) 

C  TC  =  chamber  temperature  (R) 

C  W  =  mass  in  chamber  (lb) 

C  Tdot  =  time  derivative  of  temperature  (R/s)  in  chamber 
C  Wdot  =  time  derivative  of  mass  (Ib/s)  in  chamber 
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n  n 


C  Time-dependant  parameters: 

C 

C  PO  =  gun  chamber  pressure  ^si)  -  from  XKTC  lookup  table 
C  TO  =  gun  chamber  temperature  (R)  -  from  XKTC  lookup  table 

C  Initial  Conditions: 

C 

C  PCO  =  initial  chamber  pressure  (lbJG^ft'^2) 

C  TCO  =  initial  chambCTtempCTature(R) 

WO  =  initial  mass  in  chamber  (R) 
timeO  =  start  time  (ms) 

C  timef  =  end  time  (ms) 

C- - 


IMPLICIT  DOUBLE  PREQSION  (A-H,0-Z) 

DOUBLE  PRECISION  K>I 

EXTERNAL  DDASSL 
EXTERNAL  RES30SW,DYNAM 

INTEGER  NEQ,NGJMAXORD,LIW,LRW 

PARAMETER(N  =  15) 

PARAMETER(NEQ  =3) 

PARAMETER(MAXORD  =  5) 

PARAMETEROUW  =20+NEQ) 

PARAMETER(LRW  =  40+(MAXORD+4)*NEQ+NEQ**2) 

INTEGER  INFO(N),  IDID,  IWORK(LIW),  IPAR 

DOUBLE  PRECISION  TIME,  TIMEO.TIMEF,  TIMEOUT, 

&  Y(NEQ),  YPRIME(NEQ), 

&  RTOL,  ATOL,RWORK(LRW),  RPAR 

LOGICAL  IFLAG 
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SAVEJLO 


c  Common  blocks  for  input  data. 
#include  "data.inc" 


CALLGETDATAQ 

PC0  =  PRESX(1) 

TC0  =  TEMPX(1) 

CALL  EOSW(PCO,TCO,WO) 

TIME0  =  TIMEX(1) 

TIMEF  =  TIMEX(ND  ATA) 

NPTS  =  1000 

DTIME  =  (TIMEF-TIMEO)/(NPTS) 

RTOL=  1.0D-06 
ATOL  =  O.ODO 

CALL  DYNAM(YPRIME(1),YPRIME(2),TIME0,PC0,TC0,W0) 
YPRIME(3)  =  0.0D0 


c  set  initial  conditions 


Y(1)  =  W0 
Y(2)  =  TC0 
Y(3)=PC0 

CALL  SETINFO(INFO) 

TIME  =TIME0 
TIMEOUT  =  TIME  +  DTIME 

CALL  OUT(TIME,TIMEOUT,Y,YPRIME,TCO,PCO) 
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DO  100 1  =  1,  NPTS-1 

CALL  DDASSL  (RES,  NEQ,  TIME,  Y,  YPRIME,  TIMEOUT,INFO, 
&  RTOL,  ATOL,  IDID,  RWORK,  LRW,  IWORK, 

&  LIW,  RPAR,  IPAR,  JAC) 

CALL  INTERP(TIMEX,TEMPX,PRESX,NDATA,TIME,TO,PO) 
CALL  OUT(TIME,TIMEOUT,Y,YPRIME,TO,PO) 

TIME  =TIMEOUT 
TIMEOUT  =  TIME  +  DTIME 
100  CONTINUE 

END 


SUBROUTINE  SETINFO(INFO) 

c  Pmpose:  initialize  information  array 

INTEGER  INFO(15) 

DO  1001=1, 15 
INFO(I)=0 
100  CONTINUE 

END 


c  - 

SUBROUTINE  OUT(T,TOUT,Y,YPRIME,TO,PO) 


c  Purpose:  print  out  results 

IMPLICIT  DOUBLE  PRECISION  (A-H,0-Z) 

DOUBLE  PRECISION  Y(*),YPRIME(*),PHAT 

PHAT  =  P0/Y(3) 

WRITE(6,’(1 1E20.10)')  T,Y(3)/144.0D0,Y(2),Y(1),T0,P0/144.0, 
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&  PHAT,YPRIME(1),YPRIME(2),YPRIME(3),1.0D0/PHAT 
END 


SUBROUTINE  GETDATAQ 


c  Purpose:  read  in  input  data 


CHARACTER*80  FILENAME 


#include  "data.inc" 


READ(*  *)  FILENAME 
READ(*  *)  FLOSS 
READ(*  *)  Vc 
READ(*  *)  A 
READ(*  *)  Cd 

OPEN(20,FILE=’xktc.data') 

REWIND(20) 

DO  100I=1,NMAX 

READ(20  *,END=200)  TIMEXa),  PRESK®,  TEMPXa) 


c  convert  pressure  from  psi  to  Ib^t''2 

PRESX®  =  144.0  *  PRESX® 

100  CONTINUE 

200  NDATA  =  I-1 

END 

SUBROUTINE  RES(T,Y,YPRIME, DELTA, IRES,RPARJPAR) 

DOUBLE  PREQSION  Y(*),YPRIME(*),DELTA(*) 

DOUBLE  PREQSION  P0,f  0,TC,PC,W,WDOT,PDOT,TDOT 
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EXTERNAL  EOSP,  DYNAM 


W  =Y(1) 
TC  =  Y(2) 
PC  =  Y(3) 


c  compute  dynamic  variables 


CALL  DYNAM(WDOT,TDOT,T,PC,TC,W) 


c  compute  pressure  from  equation  of  state 


CALLEOSP(PC,TC,W) 

DELTA(l)  =  YPRIME(l)  -  WDOT 
DELTA(2)  =  YPRI1V[E(2)  -  TDOT 
DELTA(3)  =  Y(3)-PC 

END 


SUBROUTINE  DYNAM(WDOT,TDOT,TIME,PC,TC,W) 

c  Describe  dynamic  equations  for  system 

DOUBLE  PRECISION  P0,T0,TC,PC,W 

DOUBLE  PRECISION  WDOT,PDOT,TDOTT>HAT,TIMEAl 


c  common  black  for  input  data 

#include  "data.inc" 

#include  "constants-inc" 

C  Obtain  gun  chamber  pressure  and  temp^ature  from  lookup  table: 
CALL  INTERP(T[MEX,PRESX,TEMPX,NDATA,TIME,PO,TO) 
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C  Compute  time  derivative  of  mass  in  chamber.  Note  that 
C  PHAT  <  0.556  implies  choked  flow.  For  P_gun  >  PC, 


IF  (PO  .GE.  PC)  THEN 

PHAT  =  PC/P0 

IF  (PHAT  .GT.  0.556D0)  THEN 

M  =  dsqrt(((-1.0D0  +  (1.0D0/PHAT)**((k-1.0D0)/k)) 
&  *2.0D0/(k-1.0D0))) 

Wdot  =  Cd  *  A  *  dsqrt(k  *  g/R)  *  (PC / dsqrt(T0))  * 
&  M  *  dsqrt(1.0D0  +  0.5D0  *  (k-l.ODO)  *  M**2) 

ELSE 

Wdot  =  Cd  *  A  *  dsqrt(  (k  *  g/R)  * 

&  ((2.0DO/(k+1.0DO))**((k+1.0DO)/(k-1.0DO)))) 

&.  *P0/dsqrt(T0) 

ENDIF 

Tdot  =  (Wdot/(W*Cv))  *  (Cp*T0*floss  -  Cv*TC) 

ELSE 

PHAT  =  P0/PC 

IF  (PHAT  .GT.  0.556DO)  THEN 

M  =  dsqrt(((-1.0D0  +  (1.0D0/PHAT)**((k-1.0D0)/k)) 
&  *2.0D0/(k-1.0D0))) 

Wdot  =  -Cd  *  A  *  dsqrt(k  *  g  /  R)  * 

&  (PO  /  dsqrt(TC))  *  M  *  dsqrt(1.0D0+0.5D0* 

&  (k-1.0D0)*M**2) 

ELSE 

Wdot  =  -Cd  *  A  *  dsqrt(  (k  *  g  /  R)* 

&  ((2.0D0/(k+1.0D0))**((k+1.0D0)/(k-1.0D0)))) 

&  *PC/dsqrt(TC) 

ENDIF 

Tdot  =  (Wdot/(W*Cv))  *  TC  *  (Cp  -  Cv) 

ENDIF 
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END 


SUBROUTINE  EOSP(PC,TC,W) 

c  equation  of  state  for  computing  chamber  pressure 

DOUBLE  PREaSION  PC,TC,W 

#include  "data.inc" 

#include  "constants.inc" 

PC  =  Ru*TC*W/(Vc*MW-b*W) 

END 

SUBROUTINE  EOSW(PC,TC,W) 

c  equation  of  state  for  computing  mass  inside  chamber 

DOUBLE  PRECISION  PC,TC,W 

#include  "data.inc" 

#include  "constants.inc" 


W  =PC*  Vc*MW/(PC*b  +  Ru*TC) 

END 

SUBROUTINE  INTERP(T^,Y,N,TINT^INT,YINT) 

- - 

c  Linear  interpolation  fiom  data  table, 
c 

c  Input: 

c  T  =  ordered  array  of  length  N  representing  the 
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c  data  domain  (must  be  either  monotonically 

c  increasing  or  descreasing). 

c  X  =  data  array  of  length  N. 

c  Y  =  data  value  of  length  N. 

c  TINT  =  interpolation  point, 
c 

c  Output: 

c  XENT  =  interpolated  X-value  at  point  TINT 
c  YINT  =  inteipolated  X-value  at  point  TINT 


DOUBLE  PRECISION  T(*),X(*),Y(*),TINT,XINT,YINT 

DOUBLE  PREaSION  THETA 

SAVEJLO 

CALL  HUNT(T,N,TINTJLO) 

IF  (JLO  .LE.  0)  THEN 
XINT  =  X(1) 

YINT  =  Y(1) 

ELSE  IF  (JLO  .GE.  N)  THEN 
XINT  =  X(N) 

YINT  =  Y(N) 

ELSE 

THETA  =  (TINT  -  T(JLO))/(T(JLO+l)-T(JLO)) 

XINT  =  X(JLO)  +  THETA  *  (X(JLO+l)  -  X(JLO)) 
YESrr  =  Y(JLO)  +  THETA  *  (Y(JLO+l)  -  Y(JLO)) 
ENDIF 

END 


SUBROUTINE  HUNT(XX,N,X,JLO) 

c  Purpose:  Search  ordered  list, 
c 

c  Source:  This  subroutine  was  copied  from  the  following  book: 
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c 

c  W.H.  Press,  B.P.  Flannery,  S.A.  Teukolsky,  W.T.  Vetterling 

c  Numerical  Recipes 

c  The  Art  of  Scientific  Computing 

c  Cambridge  University  Press,  New  York,  1986. 

c  Pages  91-92. 

c 

c  Input: 

c  XX  =  ordered  array,  must  be  eithCT  monotonically 
c  increasing  or  descreasing 

c  n  =  length  of  array  XX 
c  X  =  data  value 
c  jlo  =  initial  guess  for  index  search 
c 

c  Output 

c  jlo  =  index  value  for  xx(jlo)  <=  x  <  xx(jlo+l), 
c  jlo=0  or  jlo=N  indicates  that  x  is  out  of 

c  range. 

c - 


DOUBLE  PRECISION  XX(N),X 

LOGICAL  ASCND 

ASCND  =  XX(N)  .GT.  XX(1) 

IF  (JLO  .LE.  0  .OR.  JLO  .GT.  N)  THEN 
JLO  =  0 
JHI  =  N+1 
GOTO  3 
ENDIF 


INC  =  1 

IF  (X  .GE.  XX(JLO)  .EQV.  ASCND)  THEN 
1  JHI  =  JLO  +  INC 
IF  (JHI  .GT.  N)  THEN 
JHI  =  N+1 

ELSE  IF  (X  .GE.  XX(JHI)  .EQV.  ASCND)  THEN 
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JLO  =  JHI 
IN  =  INC  +  INC 
GOTOl 
ENDEF 
ELSE 
JHI  =  JLO 

2  JLO  =  JHI-INC 

IF  (JLO  XT.  1)  THEN 
JLO  =  0 

ELSE  IF  (X  .LT.  XX(JLO)  ,EQV.  ASCND)  THEN 
JHI  =  JLO 
INC  =  INC  +  INC 
GOTO  2 
ENDIF 
ENDIF 

3  CONTINUE 

IF  (JHI  -  JLO  .EQ.  1)  RETURN 
JM  =  (JHI  +  JLO)/2 

IF  (X  .GT.  XX(JM)  .EQV,  ASCND)  THEN 
JLO  =  JM 
ELSE 
JHI  =  JM 
ENDIF 

GOTO  3 
END 


Filename:  data.inc 


C- - 

C  Design  Parameters: 
C 
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C  Vc  =  volume  of  chamber  [ft'^3] 

C  A  =  area  of  duct  [ft^2] 

C  Cd  =  discharge  coefficient 
C  floss  =  fraction  heat  loss  from  tube 
C 

C  TIMEXO  =  table  from  XKTC  output  of  time  values  (ordered) 

C  PRESXO  =  table  from  XKTC  ou^ut  of  gun  pressure  values  (lb/ft'^2) 
C  TEMPXQ  =  table  from  XKTC  output  of  gun  temperature  values  (R) 
C- - 


PARAMETER(NMAX  =5000) 

DOUBLE  PRECISION  TIMEX(NMAX),  PRESX(NMAX),  TEMPX(NMAX) 
DOUBLE  PRECISION  FLOSS,A,Vc,Cd 
INTEGER  NDATA 

COMMON  /XKTCl/  TIMEX,  PRESX,  TEMPX 
COMMON  /XKTC2/  NDATA 
COMMON  /PARAM/  FLOSS,  Vc,A,Cd 


C- - 

C  Original  problem: 

C 

C  floss  =  0.85D0 
C  Vc  =0.005787D0 
C  A  =0.000608D0 
C  Cd  =0.62D0 
C- - 


Filename:  constants.inc 


C- - 

C  Constants: 

C 

C  R  =  gas  constant  [ft-lb/(lb  R)] 

C  MW  =  molecular  weight  [Ib/lb-mole] 
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a 


C  g  =  gravitational  acceleration  [ft/sec'^2] 

=  constant  volume  heat  capacity  [ft-lb/(lb  R)] 

=  constant  pressure  heat  capacity  [ft-lb/(lb  R)] 

C  J  =  convertion  factor  [ft-lb/Btu] 

C  k  =  ratio  of  Cp  to  Cv 
C  Ru  =  universal  gas  constant  [ft-lb/(lb  R  mol)] 

C  b  =  specific  covolume  [ft^3/lb] 

C- - 

DOUBLE  PRECISION  R,MW,g,J,k,Cp,Cv,Ru,b 

PARAMETER(R  =  66.09D0) 
PARAMETER(MW  =  23.36D0) 

PARAMETER(g  =  32.2D0) 

PARAMETER(J  =778.0D0) 

PARAMETER(k  =  1.243D0) 

PARAMETER(Cp  =  336.7D0) 

PARAMETER(Cv  =  27O.698D0) 
PARAMETER(Ru  =  1544.0D0) 

PARAMETER(b  =O.01649D0) 
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2  ADMINISTRATOR 
ATTN  DUC  DDA 

DEFENSE  TECHNICAL  INFO  CTR 
CAMERON  STATION 
ALEXANDRIA  VA  22304-6145 

1  DIRECTOR 

ATTN  AMSRL  OP  SD  TA 
US  ARMY  RESEARCH  LAB 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

3  DIRECTOR 

ATTN  AMSRL  OP  SD  TL 
US  ARMY  RESEARCH  LAB 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

1  DIRECTOR 

ATTN  AMSRL  OP  SD  TP 
US  ARMY  RESEARCH  LAB 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 


ABERDEEN  PROVING  GROUND 

2  DIR  USARL 

ATTN  AMSRL  OP  AP  L  (305) 
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1  HQDA 

ATTN  SARD  TR  MS  K  KOMINOS 
PENTAGON 

WASHINGTON  DC  20310-0103 
1  HQDA 

ATTN  SARD  TR  DR  R  CHATT 
PENTAGON 

WASHINGTON  DC  20310-0103 

1  CHAIRMAN 

DOD  EXPLOSIVES  SAFETY  BD 
HOFFMAN  BLDG  1  RM  856  C 
2461  EISENHOWER  AVE 
ALEXANDRIA  VA  22331-0600 

1  HQS  US  ARMY  MATERIEL  CMD 

ATTN  AMCICP  AD  M  FISETTE 
5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

1  US  ARMY  BMDS  CMD 

ADVANCED  TECHLGY  CTR 
PO  BOX  1500 

HUNTSVILLE  AL  35807-3801 

1  OFC  OF  THE  PRODUCT  MGR 

155MM  HOWITZER  M109A6  PALADIN 
ATTN  SFAE  AR  HIP  IP 
MR  R  DE  KLEINE 
PICAHNNY  ARSENAL  NJ 
07806-5000 

3  PMAFAS 

ATTN  SFAE  ASM  AF  E 
LTC  A  ELLIS 
TKURIATA 
J  SHIELDS 

PICAHNNY  ARSENAL  NJ 
07806-5000 

1  PMAFAS 

ATTN  SFAE  ASM  AF  Q  W  WARREN 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  CDR  US  ARMY  ARDEC 

PROD  BASE  MDRNZTN  AGENCY 
ATTN  AMSMC  PBM  A  SIKLOSI 
PICATINNY  ARSENAL  NJ 
07806-5000 


1  CDR  US  ARMY  ARDEC 

PROD  BASE  MDRNZTN  AGENCY 
ATTN  AMSMC  PBM  E  L  LAIBSON 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  PM  PEO  ARMAMENTS 

TANK  MAIN  ARMAMENT  SYSTEM 
ATTN  AMCPM  TMA 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  PM  PEO  ARMAMENTS 

TANK  MAIN  ARMAMENT  SYSTEM 
ATTN  AMCPM  TMA  105 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  PM  PEO  ARMAMENTS 

TANK  MAIN  ARMAMENT  SYSTEM 
ATTN  AMCPM  TMA  120 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  PM  PEO  ARMAMENTS 

TANK  MAIN  ARMAMENT  SYSTEM 
ATTN  AMCPM  TMA  AS  H  YUEN 
PICATINNY  ARSENAL  NJ 
07806-5000 

2  CDR  US  ARMY  ARDEC 
ATTN  SMCAR  CCH  V 
CMANDALA 

E FENNELL 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  CDR  US  ARMY  ARDEC 
ATTN  SMCAR  CCH  T 
L  ROSENDORF 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  CDR  US  ARMY  ARDEC 
ATTN  SMCAR  CCS 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  CDR  US  ARMY  ARDEC 

ATTN  SMCAR  AEE  J  LANNON 
PICATINNY  ARSENAL  NJ 
07806-5000 
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COMMANDER 

US  ARMY  ARDEC 

ATTN  SMCAR  AEE  WW 

MMEZGER 

J  PINTO 

DWIEGAND 

PLU 

CHU 

PICATINNY  ARSENAL  NJ 
07806-5000 

CDR  US  ARMY  ARDEC 
ATTN  SMCAR  AEE  B 
D  DOWNS 
S  EINSTEIN 
S  WESTLEY 
S  BERNSTEIN 
J  RUTKOWSKI 
B  BRODMAN 
P  O’REILLY 
RCmiNaONE 
PHUI 

JO’REELLY 

PICATINNY  ARSENAL  NJ 
07806-5000 

COMMANDER 
US  ARMY  ARDEC 
ATTN  SMCAR  AES 
SKAPLOWnZ 
PICATINNY  ARSENAL  NJ 
07806-5000 

COMMANDER 
US  ARMY  ARDEC 
ATTN  SMCAR  FSA  T 
M  SALSBURY 
PICATINNY  ARSENAL  NJ 
07806-5000 

COMMANDER 
US  ARMY  ARDEC 
ATTN  SMCAR  HFM 
EBARRIERES 
PICATINNY  ARSENAL  NJ 
07806-5000 

COMMANDER 

US  ARMY  ARDEC 

ATIN  SMCAR  FSC  G  FERDINAND 

PICATINNY  ARSENAL  NJ 

07806-5000 


1  COMMANDER 

US  ARMY  ARDEC 
ATTN  SMCAR  FS  T  GORA 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
ATTN  SMCAR  FS  DH  J  FENECK 
PICATINNY  ARSENAL  NJ 
07806-5000 

3  COMMANDER 

US  ARMY  ARDEC 
ATTN  SMCAR  FSS  A 
RKOPMANN 
B  MACHEK 
L  BINDER 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 

ATTN  SMCAR  FSN  N  K  CHUNG 
PICATINNY  ARSENAL  NJ 
07806-5000 

2  DIR  BENET  LABS 
ATTN  SMCAR  CCB  RA 
GP  O’HARA 

G  A  PFLEGL 

WATERVLIET  NY  12189-4050 

1  DIR  BENET  LABS 

ATTN  AMSTA  AR  CCB  T 
SSOPOK 

WATERVLIET  NY  12189-4050 

1  DIR  BENET  LABS 
ATTN  SMCAR  CCB  S 
FHEISER 

WATERVLIET  NY  12189-4050 

2  CDR  US  ARMY  RSRCH  OFC 
ATTN  TECHNICAL  LIBRARY 
DMANN 

PO  BOX  12211 

RESEARCH  TRIANGLE  PARK  NC 
27709-2211 
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1  CDR  USACECOM 

R&D  TECHNICAL  LIBRARY 
ATTN  ASQNC  ELC  IS  L  R 
MYERCTR 

FT  MONMOUTH  NJ  07703-5301 
1  CMDT 

US  ARMY  AVIATION  SCHOOL 
ATTN  AVIATION  AGHSfCY 
FT  RUCKER  AL  36360 

1  PM 

US  TACOM 

ATTN  AMCPM  ABMS  T  DEAN 
WARREN  MI  48092-2498 

1  PM  US  TANK  AUTOMOTIVE  CMD 
nGHlTNG  VEHICLE  SYSTEMS 
ATTN  SFAE  ASM  BV 
WARREN  MI  48397-5000 

1  PM  ABRAMS  TANK  SYSTEM 
ATTN  SFAE  ASM  AB 
WARREN  MI  48397-5000 

1  DIR  HQ  TRAC  RPD 
ATTNATCDMA 
FT  MONROE  VA  23651-5143 

1  COMMANDANT 

US  ARMY  CMD  &  GEN 
STAFF  COLLEGE 
FT  LEAVENWORTH  KS  66027 

1  COMMANDANT 

US  ARMY  SPCL  WARFARE  SCHL 
ATTN  REV  AND  TRNG  LIT  DIV 
FT  BRAGG  NC  28307 

1  COMMANDER 

RADFORD  ARMY  AMMO  PLANT 
ATTN  SMCAR  QA  HI  LIB 
RADFORD  VA  24141-0298 

1  COMMANDER 
US  ARMY  NGIC 
ATTN  AMXST  MC  3 
220  SEVENTH  STREET  NE 
CHRLTTESVLLE  VA  22901-5396 


2  COMMANDANT 
US  ARMY  FAC&S 
ATTN  ATSF  CO  MW  E  DUBUSKY 
ATSF  CN  P  GROSS 
FTSDXOK  73503-5600 

1  COMMANDANT 

US  ARMY  FIELD  ARTUIY 
CTR  &  SCHOOL 
ATTN  ATSF  CN 
FT  SILL  OK  73503-5600 

1  CMDT  US  ARMY  ARMOR  SCHOOL 
ARMOR  AGENCY 

ATTN  ATZK  CD  MS 
M  FALKOVrrCH 
FT  KNOX  KY  40121-5215 

2  CDR  NAVAL  SEA  SYSTEMS  CMD 
ATTN  SEA  62R 

SEA  64 

WASH  DC  20362-5101 

1  CDR 

US  ARMYARDEC 
ATTN  SMCAR  FSA  F 
S  FLOROFF 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  CDR  NAVAL  AIR  SYSTEMS  CMD 

ATTN  AIR  954  TECH  LIB 
WASHINGTON  DC  20360 

4  CDR  NAVAL  RSRCH  LAB 

ATTN  TECHNICAL  LIBRARY 
CODE  4410 
K  KAILASANATE 
J  BORIS 
EORAN 

WASHINGTON  DC  20375-5000 

1  OFFICE  OF  NAVAL  TECHLGY 

ATTN  ONT  213  D  SIEGEL 
800  N  QUINCY  ST 
ARLINGTON  VA  22217-5000 
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7  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
ATTN  T  C  SMITH 
KRICE 
S  MITCHELL 

5  PETERS 

J  CONSAGA 
CGOTZMER 
TECHNICAL  LIBRARY 
INDIAN  HEAD  MD  20640-5000 

1  COMMANDER 

NAVAL  SURFACE  WARFARE  CIR 
ATTN  CODE  730 
SILVER  SPRING  MD 
20903-5000 

1  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
ATTN  CODE  R  13  R  BERNECKER 
SILVER  SPRING  MD 
20903-5000 

1  DIRECTOR 

US  ARMY  TRAC  FT  LEE 
ATTN  ATRC  L  MR  CAMERON 
FORT  LEE  VA  23801-6140 

1  COMMANDER 

US  ARMY  BELVOm  RD&E  CTR 
ATTN  STRBE  WC 
FT  BELVOIR  VA  22060-5606 

1  US  ARMY  RSCH  DEVELOPMENT 

6  STNDRDZTN  GROUP  UK 
ATTN  DR  ROY  E  RICHENBACH 
PSC  802  BOX  15 

FPO  AE  09499-1500 

1  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
ATTN  CODE  G30  GUNS 
&  MUNITIONS  DIV 
DAHLGRENVA  22448-5000 

1  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
ATTN  CODE  G32  GUNS 
SYSTEMS  DIV 
DAHLGRENVA  22448-5000 


1  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
ATTN  CODE  G33  T  DORAN 
DAHLGRENVA  22448-5000 

1  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
ATTN  CODE  E23 
TECHNICAL  LIBRARY 
DAHLGRENVA  22448-5000 

2  COMMANDER 

NAVAL  AIR  WARFARE  CTR 
ATTN  CODE  388 
CF  PRICE 
T BOGGS 

CHINA  LAKE  CA  93555-6001 

2  COMMANDER 

NAVAL  AIR  WARFARE  CTR 
ATTN  CODE  3895 
TPARR 
RDERR 

CHINA  LAKE  CA  93555-6001 

1  COMMANDER 

NAVAL  AIR  WARFARE  CTR 
INFORMATION  SCIENCE  DIV 
CHINA  LAKE  CA  93555-6001 

1  COMMANDING  OFFICER 

NAVAL  UNDERWATER  SYSTEMS  CTR 
ATTN  CODE  5B331  TECH  LIBRARY 
NEWPORT  RI  02840 

1  AFOSR  NA 

ATTN  J  TISHKOFF 
BOLLING  AFB  DC 
20332-6448 

1  OLAC  PL  TSTL 

ATTN  D  SHIPLETT 
EDWARDS  AFB  CA  93523-5000 

3  ALLSCF 
ATTN  J  LEVINE 
L  QUINN 

T  EDWARDS 

EDWARDS  AFB  CA  93523-5000 
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WLMNAA 
ATTN  B  SIMPSON 
EGLINAFBFL  32542-5434 

WL  MNME 

ENERGETIC  MATERIALS  BR 
2306  PERIMETER  RD 
STE  9 

EGLINAFBFL  32542-5910 


DIR  SANDIA  NATL  LABS 
ATTN  R  CARLING 
COMBUSTION  RSRCH  FACILITY 
LIVERMORE  CA  94551-0469 

DIR  SANDIA  NATL  LABS 
ATTN  8741  G  A  BENEDITTI 
PO  BOX  969 

LIVERMORE  CA  94551-0969 


WLMNSH 

ATTN  R  DRABCZUK 

EGLINAFBFL  32542-5434 

NASA  LANGLEY  RSRCH  CTR 
ATTN  M  S  408 
W  SCALLION 
DWTTCOFSKI 
HAMPTON  VA  23605 

CENTRAL  INTLGNCE  AGENCY 
OFC  OF  THE  CNTRL  RFRNCES 
DISSEMINATION  BRANCH 
ROOM  GE  47  HQS 
WASHINGTON  DC  20502 

CENTRAL  INTLGNCE  AGENCY 
ATTN  J  BACKOFEN 
NHB  ROOM  5N01 
WASHINGTON  DC  20505 

SDIOTNI 

ATTN  L  H  CAVENY 
PENTAGON 

WASHINGTON  DC  20301-7100 

SDIO  DA 
ATTN  E  GERRY 
PENTAGON 

WASHINGTON  DC  21301-7100 


2  DIRLLNL 
ATTN  L  355 
A  BUCKINGHAM 
M  FINGER 
PO  BOX  808 

LIVERMORE  CA  94550-0622 
1  DIR 

LOS  ALAMOS  NATL  LAB 
ATTN  T3  D  BUTLER 
PO  BOX  1663 
LOS  ALAMOS  NM  87544 

1  DIR  LOS  ALAMOS  NATL  LAB 

ATTN  M  DIVISION 

B  CRAIG 

PO  BOX  1663 

LOS  ALAMOS  NM  87544 

2  BATTELLE 
ATTN  TWSTIAC 
V  LEVIN 

505  KING  AVE 
COLUMBUS  OH  43201-2693 

1  BATTELLE  PNL 

ATTN  M  C  C  BAMPTON 
PO  BOX  999 
RICHLAND  WA  99352 


HQDNA 
ATTN  D  LEWIS 
A  FAHEY 

6801  TELEGRAPH  RD 
ALEXANDRIA  VA  22310-3398 


THE  UNIV  OF  AUSTIN  TEXAS 
INSTITUTE  FOR  ADVANCED  TECHLGY 
ATTNTMKIEHNE 
4030  2  W  BRAKER  LANE 
AUSTIN  TX  78759-5329 


DIR  SANDIA  NATL  LABS 
ATTN  M  BAER 
DEPARTMENT  1512 
PO  BOX  5800 

ALBUQUERQUE  NM  87185 


1  INST  OF  GAS  TCHNLGY 
ATTN  D  GIDASPOW 
3424  S  STATE  ST 
CHICAGO  IL  60616-3896 
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2  CPIAJHU 

ATTN  H  J  HOFFMAN 
T  CHRISTIAN 

10630  LTITLE  PATUXENT  PKWY 
SUITE  202 

COLUMBIA  MD  21044-3200 

1  CALFORNIA  INST  OF  TECH 

JET  PROPULSION  LABORATORY 
ATTN  L  STRAND  MS  125  224 
4800  OAK  GROVE  DRIVE 
PASADENA  CA  91109 

3  GEORGIA  INST  OF  TCHNLGY 
SCHL  OF  AEROSPACE  ENGNRNG 
ATTN  E  PRICE 
WCSTRAHLE 

BTZINN 

ATLANTA  GA  30332 

1  UNIVERSITY  OF  MARYLAND 

ATTN  DR  J  D  ANDERSON 
COLLEGE  PARK  MD  20740 

1  MASSACHUSETTS  INST  OF 

TECHNOLOGY 
ATTN  T  TOONG 
77  MASSACHUSETTS  AVE 
CAMBRIDGE  MA  02139-4307 

1  UNIVERSITY  OF  ILLINOIS 

DEPT  OF  MECH  ENG 
ATTNHKRIER 
144MEB  1206  W  GREEN  ST 
URBANA IL  61801 

1  UNIVERSITY  OF  MINNESOTA 

DEPT  OF  MECHNCL  ENGNRNG 
ATTN  E  FLETCHER 
MINNEAPOUS  MN  55455 

3  PENNSYLVANIA  STATE  UNIV 

DEPT  OF  MECHNCL  ENGNRNG 
ATTN  K  KUO 
CMERKLE 
STHYNELL 
VYANG 

UNIVERSITY  PARK  PA  16802 


1  AFELM  THE  RAND  CORP 

ATTN  LIBRARY  D 
1700  MAIN  STREET 
SANTA  MONICA  CA  90401-3297 

1  ARROW  TECHLGY  ASSOC  INC 
ATTN  W  HATHAWAY 

PO  BOX  4218 

SOUTH  BURLINGTON  VT  05401-0042 

2  AAI  CORPORATION 
ATTN  JFRANKLE 
D  CLEVELAND 

PO  BOX  126 

HUNT  VALLEY  MD  21030-0126 

8  ALLIANT  TECHSYSTEMS  INC 
ATTN  R  E  TOMPKINS 
J KENNEDY 
JBODE 
CCANDLAND 
L  OSGOOD 
RBURETTA 
R BECKER 
M  SWENSON 
600  SECOND  ST  NE 
HOPKINS  MN  55343 

1  ELI  FREEDMAN  AND  ASSOCIATES 
ATTN  E  FREEDMAN 
2411  DIANA  RD 
BALTIMORE  MD  21209-1525 

1  GENERAL  APPLIED  SCIENCES  LAB 

ATTN  J  ERDOS 
77  RAYNOR  AVE 
RONKONKAMANY  11779-6649 

1  GENERAL  ELECTRIC  COMPANY 

TACTICAL  SYSTEM  DEPT 
ATTN  JMANDZY 
100  PLASTICS  AVE 
PITTSFIELD  MA  01201-3698 

1  HTRl 

ATTN  M  J  KLEIN 
10  W  35TH  STREET 
CHICAGO  IL  60616-3799 
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4  HERCULES  INC 
ATTNLGIZZr 
D  A  WORRELL 
WJ  WORRELL 
C CHANDLER 

RADFORD  ARMY  AMMO  PLANT 
RADFORD  VA  24141-0299 

2  HERCULES  INC 

ATTN  WILLIAM  B  WALKUP 
THOMAS  F  FARABAUGH 
ALLEGHENY  BALLISTICS  LAB 
PO  BOX  210 

ROCKET  CENTER  WV  26726 

1  HERCULES  INC 

ATTN  R  CARTWRIGHT 
AEROSPACE 
100  HOWARD  BLVD 
KENVILLENJ  07847 

1  HERCULES  INC 

ATTN  B  M  RIGGLEMAN 
HERCULES  PLAZA 
WILMINGTON  DE  19894 

1  MARTIN  MARIETTA  DEFENSE  SYSTEMS 

ATTN  GEORGE  KEELER 
RM2260 

100  PLASTICS  AVE 
PITTSFIELD  MA  01201 

1  MBR  RESEARCH  INC 

ATTN  DR  MOSHE  BEN  REUVEN 
601  EWING  ST  SUHE  C  22 
PRINCETON  NJ  08540 

1  OLIN  CORPORATION 

ATTN  F  E  WOLF 
BADGER  ARMY  AMMO  PLANT 
BARABOOWI  53913 

3  OLIN  ORDNANCE 

ATTN  E  J  KIRSCHKE 
A  F  GONZALEZ 
D  W  WORTHINGTON 
PO  BOX  222 

ST  MARKS  FL  32355-0222 


1  OLIN  ORDNANCE 

ATTN  H  A  MCELROY 
10101  9TH  STREET  NORTH 
ST  PETERSBURG  FL  33716 

1  PAUL  GOUGH  ASSOC  INC 

ATTN  P  S  GOUGH 

1048  SOUTH  ST 
PORTSMOUTH  NH  03801-5423 

1  PHYSICS  INTERNATIONAL  LIBRARY 
ATIN  H  WAYNE  WAMPLER 

PO  BOX  5010 

SAN  LEANDRO  CA  94577-0599 

2  PRINCETON  CMBSTN  RSRCH  LABS  INC 
ATTN  N  A  MESSINA 

NMER 

PRINCETON  CORPORATE  PLAZA 
11  DEERPARK  DR  BLDG  IV  SUITE  119 
MONMOUTH  JUNCTION  NJ  08852 

1  TEXTRON  DEFENSE  SYSTEM 
ATTN  A  PATRICK 

2385  REVERE  BEACH  PKWY 
EVERETT  MA  02149-5900 

3  ROCKWELL  INTRNTNL 
ATTN  BA08 

J  FLANAGAN 

JGRAY 

R  B  EDELMAN 

ROCKETDYNE  DIV 

6633  CANOGA  AVE 

CANOGAPARKCA  91303-2703 

2  ROCKWELL  INTRNTNL  SCIENCE  CTR 
ATTN  DR  S  CHAKRAVARTHY 

DR  S  PALANISWAMY 

1049  CAMINO  DOS  RIOS 
THOUSAND  OAKS  CA  91360 

1  SAIC 

ATTN  M  PALMER 
2109  AIR  PARK  RD 
ALBUQUERQUE  NM  87106 

1  SOUTHWEST  RSRCH  INSTITUTE 

ATTN  J  P  RBEGEL 
6220  CULEBRA  ROAD 
SAN  ANTONIO  TX  78228-0510 
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1  SVERDRUP  TECHLGY  INC 

ATTN  DR  JOHN  DEUR 
2001  AEROSPACE  PARKWAY 
BROOK  PARK  OH  44142 

3  THIOKOL  CORPORATION 

ATTNRWILLER 
R  BIDDLE 
TECH  LIBRARY 
ELKTON  DIVISION 
PO  BOX  241 

ELKTON  MD  21921-0241 

1  VERTTAY  TECHLGY  INC 

ATTN  E  FISHER 
4845  MILLERSPORT  HWY 
EAST  AMHERST  NY  14501-0305 

1  UNIVERSAL  PROPULSION  COMPANY 

ATTN  H  J  MCSPADDEN 
25401  NORTH  CENTRAL  AVE 
PHOENIX  AZ  85027-7837 

1  SRI  INTERNATIONAL 

ATTN  TECH  LIBRARY 
PROPULSION  SCIENCES  DIV 
333  RAVENWOOD  AVE 
MENLO  PARK  CA  94025-3493 


ABERDEEN  PROVING  GROUND 

1  CDR  USA  ATC 
ATTN  SIECS  LI 
R  HENDRICKSEN 
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NO.  OF 
COPIES 


54 


ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


ABERDEEN  PROVING  GROUND 


DIRUSARL 

ATTN  AMSRL  WT  P  A  HORST 
AMSRLWTPA 
T  MINOR 
ABIRK 
K  WHITE 

L-M  CHANG  (5  CP) 

J COLBURN 
P  CONROY 
G KELLER 
DKOOKER 
MNUSCA 
TROSENBERGER 
A  BRANT 
S  HOWARD  (5  CP) 

D  KRUCZYNSKI 
CRUTH 
ISTOBIE 
A  WILLIAMS 
AMSRL  WT  PB 
P  PLOSTINS 
AMSRL  WT  PC 
G  ADAMS 
R  ANDERSON 
MMILLER 
AMSRL  WTT 
W  MORRISON 
AMSRL  WTTB 
RFREY 
AMSRL  WTTC 
W  DE  ROSSET 
F GRACE 
B  SORENSEN 
AMSRL  WT  ID 
A  DIETRICH 
KFRANK 
AMSRL  WTWA 
B  MOORE 
H  ROGERS 
ABARAN 
AMSRL  WTWB 
F BRANDON 


AMSRL  SC 
W  MERMAGEN 
AMSRL  see 
WSTUREK 
AMSRL  SC  CC 
ACELMINS 
J  GROSH  (5  CP) 
AMSRL  SC  I 
M  HIRSCHBERG 
RKASTE 
AMSRL  SC  S 
AMARK 
AMSRL  SL 
J  SMITH 
AMSRL  SL  I 
D  HASKELL 
AMSRL  SL  BG 
DKIRK 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  qualify  of  the  reports  it  publishes.  Your  comments/answers 
to  the  itemsi^questions  below  will  aid  us  in  om:  efforts. 

1.  ARL  Report  Numbei/Author  ARL-TR-1094  (Howard) _ Date  of  Report  June  1996 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  piupose,  related  project,  or  other  area  of  interest  for  which  the  report 

will  be  used.) _ 


4.  Specifically,  how  is  tiie  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.) _ 


Organization 

CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 

City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the 
Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


